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Abstract Suppose that a large number of parameterized trajectories 7 of a dynamical system evolving in R are 
stored in a database. Let ij C R w denote a parameterized path in Euclidean space, and let || • || denote a norm on 
the space of paths. In this paper, we define data structures and indices for trajectories and give algorithms to answer 

queries of the following forms: , lf 

Query 1. Given a path 77, determine whether 77 occurs as a subtrajectory of any trajectory 7 from the database. 

so, return the trajectory; otherwise, return null. 

Query 2. Given a path 77, return the trajectory 7 from the database which minimizes the norm 

11*1 — ill- 


1 Queries about trajectories 

Suppose that a large number of parameterized trajectories 7 of a dynamical system evolving in R are 
stored in a database. Let 77 C R W denote a parameterized path in Euclidean space, and let || • || denote a 
norm on the space of paths to be specified later. In this paper, we define a data structure and indices to 
represent trajectories of dynamical systems and sketch algorithms to answer queries of the following forms: 
Query 1. Given a path 77, determine whether 77 occurs as a subtrajectory of any trajectory 7 from the 
database. If so, return the trajectory; otherwise, return null. 

Query 2. Given a path 77, return the trajectory 7 from the database which minimizes the norm 

11*7 — Til- 

The paper is a successor to [2], which describes the data structure to store trajectories which is used here. 

Efficient algorithms to answer these type of queries should prove useful for a number of applications. 
As an example, consider the path-planning problem for a robotic arm. Suppose that a large number of 
feasible trajectories of the robotic arm have been stored in a database. Let 77 be the desired path of the 
arm. It is not necessary that 7, itself be a feasible trajectory. Query 2 would return the feasible trajectory 
7 of the arm which is closest to the desired path i). 

As another example, consider a database containing control trajectories for an aircraft. Assume that 
those trajectories which enter into an unstable control regime somewhere along their flight path are tagged. 
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Let rj denote a measured portion of the flight path. Then Query 1 would return the nearest full control 
trajectory in the database, which includes information about the stability of the trajectory. More generally, 
one could imagine retrieving from the database those stable trajectories which avoid a given obstacle, such 
as a turbulent region of space. In other words, the query could be used as part of a supervisory control 
system and be viewed as a means of extracting qualitative or summary information about the control 
system. 

The data structure we use to represent trajectories is closely related to hashing methods for curves that 
have been used in computer vision; see [7] and [8]. A related means of extracting qualitative information 
from dynamical systems is described in [1]. We are concerned in this paper with data structures and 
indexing for object oriented databases consisting of trajectories. For general methods of indexing in object 
oriented databases see [3], [4], [5], and [6]. 

In Section 2 we review the relevant facts about trajectories of differential equations and define different 
data structures to store trajectories. In Section 3, we show how these data structures can be used to answer 
the queries above. Section 4 contains some concluding remarks. 

2 Paths, trees and vector fields 

In this section, we describe a data structure for paths following [2]. The point of view is to assume that 
the path arises from a trajectory of a differential equation and to base the data structure upon the initial 
value problem for the differential equation. 

We begin by recalling some basic facts and definitions about trajectories of differential equations. Let 
Dfx = dfdxp. A vector field 

N 

E='£a»D„ 

M— 1 

on R* is determined by specifying N functions 

a„ : R n — * R 

We also denote the vector field by E a . A parameterized path 

7 : [^V 1 ] C R — ► K N 

is called a trajectory of the dynamical system 


i(t) = E a (x(t)) (1) 

in case it is the unique solution of the initial value problem 

m = E a (z(t)), *(<°) = 7 (t°). (2) 

We define the vector field/reference point representation or VEFREP of a path r/ to be the pair (E } R), 
consisting of a vector field E and a reference or initial point R ) where the trajectory is the solution of the 
initial value problem 

x(t) = E(x(t)), x(t°) = R. 

Note that this representation is not unique. Indeed, several different vector fields could have a given spatial 
curve as a trajectory, while any point along the spatial curve could serve as the initial value. 

We now give an algorithm whose input is a parameterized path 

*?:[*V]C R— *R n , 
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and whose output is a labeled, rooted binary tree. We assume for convenience that t° — 0 and t — 1 , 
if not, we can reparameterize. We do not assume that r] is a trajectory of the dynamical system ( 1 ). To 
define the tree, we first fix a tolerance e > 0 . The tree we define is a subset of the complete rooted binary 
tree. There are 2 k children at height k from the root: number them left to right from 1 to 2 fc . We assign 
two labels to the jth node v from the left at height k: 

K (v) = (i/2 k ) (n(^) - )) € rW 

and > . | \ 

>W = n( 1 f t -)en" 

We use the following stopping criterion to grow the tree. If a node has children v and v with labels k and 
k', respectively, and if \\k - k ' || < e, then the nodes u and v f are leaves. Here || • || denotes the Euclidean 
norm. We denote by T(t]) the tree that arises in this fashion. This tree has a simple interpretation, the 
6 labels represent points on the path 77, while the k labels represent approximate tangent vectors at those 
points. The tree is grown until the difference between two adjacent tangent vectors is uniformly small. 

Using the tree T(r]) ) we now define a vector field E(rf). The vector field E{r\) is simply the vector field 
which interpolates the labels (0(u), 77(1;)) 

E(7]){$(v)) = k(v), (3) 

for all leaves v in T(r}). Recall that 0 (v) is the point on the curve 77 corresponding to the node u, and k(u) 
is the approximate tangent to the curve at that point. 

For some applications, it is better to impose an upper bound on the degree of the interpolating functions. 
Let q denote this bound. In this case, we can define the vector field E{r\) by requiring that the coefficients 
bj minimize the quantity 

£ \\E{tiMv)) - k(v)\\, ( 4 ) 

leaves v 

where the minimum is over vector fields with interpolating functions of degree less than or equal to q. 

We conclude this section by defining a specific point 7), associated with a parametrized path 

rj : [f 0 ,* 1 ] C R — K N 

Let T(tj) the corresponding tree and E(tj) the associated vector field. Consider the trajectory defined by 
the initial value problem 

X(t) = E a (x(t)), *( 0 ) = 77 (f°)- 

Let 7 denote this trajectory. In general, 7 is only an approximation to the path rj. Define R(r}) € IL as 
follows: if the path 7 and the unit sphere in R/^ intersect precisely once, let R(f]) denote this intersection, if 
they intersect several times, let /?( 77) denote the intersection which occurs first when the list of intersections 
is ordered in lexigraphical order; otherwise, let /?( 77) denote the closest point between the unit sphere and 
the trajectory 7. 

3 Query Algorithms 

Let rj C denote a path. In this section, we define an index I(t)) that can be used for storing and 
accessing the path 7. First, fix injective functions 


hn :R n 


{ 1 , 2 ,...}, 
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for each n ~ 1 , 2 ,... Given a path 17, first compute its VEFREP (£(77), R(tj)) from its rectifying tree T(rj ). 
Assume the tree T(rj) has K leaves. Next, view the coefficients of the vector field £(77) as a K • N vector, 
so that the pair (£(77), R(t})) has ( K + 1 )N components. Then the hash index H(rj) associated with 77 is 
defined by 

H ( T l) - h(K+i)N(E(r)),R(ri)). 

We can now assign indices I(r}) to trajectories sequentially: use the hash index 7/(77) to determine 
whether the path 77 has an index assigned to it. If so, use that index; if not, use the next available index. 
Suppose that 71, ... ,7 p are trajectories of the dynamical system 

i(f) = £' a (x(f)), 

as a ranges over some parameter space. To each such trajectory 7, let 

(E( 7), *(7)) 

denote its VEFREP representation. Given a parameterized path 77, Algorithm 1 below returns the trajec- 
tory 7 from the database which contains a segment equal to the path 17. If there is no such trajectory, null 
is returned. 

Algorithm 1. The input is a parameterized path 77, and the output is the trajectory 7 from the database 
answering Query 1 . Fix e > 0 and q > 1 . 

Step 1 . This step is a precomputation. For each trajectory 7*, i — compute its VEFREP 

representation (£(7^), #(7,-)). This depends upon q and e. 

Step 2. Given a query path 77, compute its rectifying tree T(? 7). This depends upon e. Using T(r}) and 
Equation 4 , compute its VEFREP representation (£(77), #(77)). This depends upon q. 

Step 3 . Using the VEFREP (£(77), 72(77)), compute the hash index 7/(77). If there is an index in the row 
£(77) the index table, retrieve the VEFREP representation (£, 72 ) corresponding to this index; 
otherwise, return null. 

Step 4 . If Step 3 yielded a VEFREP (£, 72 ), return the trajectory 7 which is the solution to the initial 
value problem 

i(t) = £(*(<))> *( 0 ) = 72 ; 

otherwise, return null. 

Theorem 3.1 Assume that the database contains n trajectories. Algorithm 1 answers Query 1 in time 

0 ( 1 ). 

Easy modifications of Algorithm 1 can be used to answer Query 2 in time 0 (n). 

4 Conclusion 

In this paper, we have described preliminary work concerned with queries of databases containing trajec- 
tories of differential equations. Trajectories of differential equations have many different representations. 
For the types of queries considered here, we have chosen to represent parameterized trajectories 

7 : [f 0 .* 1 ] C R — + R N 


3 Q 



by a pair, consist.ng of a vector field E on R* with polynomial coefficients and a point R € R N such that 
the trajectory is the solution of the initial value problem: 


i(t) = £(*(*)), x(i°) = R. 


We call this a VEFREP representation. Using the VEFREP representation, we 
/( 7 ) and algorithms to answer queries which retrieve subtrajectories and close 


have introduced an index 
by trajectories of a given 


query trajectory. 


References 


[1] H. Abelson and G. J. Sussman, Dynamicists’ Workbench I: 
Experiments,” R. Grossman (editor), Symbolic Computation: 

SIAM, 1989. 


“Automatic Preparation of Numerical 
Applications to Scientific Computing , 


[2] R. Grossman, “Querying Databases of Trajectories of Differential Equations I: Data Structures or 
Trajectories,” to appear in Proceedings of the 23rd Hawaii International Conference on Systems Sci- 
ences, IEEE, 1990. 


[3] W. Kim, J. Banerjee, H. T. Chou, J. F. G. Garze, and D. Woelk, “Composite Object Support in an 
Object Oriented Database System,” Proceedings of OOPSLA 198 7. 


[4] W. Kim, K-C Kim, and A. Dale, 
Oriented Concepts, Databases , and 
York, 1989. 


Indexing Techniques for Object-Oriented Databases,” in Obgect- 
Applications, W. Kim and F. H. Lochovsky, editors, ACM, New 


[5] D. Maier and J. Stein, “Indexing in an Object-Oriented DBMS,” Proceedings of the 1986 International 
Workshop on Object-Oriented Database Systems, Pacific Grove, California, 1986. 


[6] N. Paton and P. M. D. Gray, “Identification of Database Objects by Key, in Advances in Object 
Oriented Database Systems, K. R. Dittrick, editor, Springer- Verlag, Berlin, 1988, pp. 280-285. 


[71 J. T. Schwartz and M. Sharir, “Identification of Partially Obscured Objects in Twc > Dimensions by 
Matching of Noisy Characteristic Curves” International J. Robotics Research, 6 (1987), ZV 

[8] H. Wolfson, “On Curve Matching,” Proceedings of Workshop on Computer Vision, Miami Beach, 
IEEE, 1987, 307-310. 


39 


